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Abstract. - Coherent macroscopic tunneling of a Bose-Einstein condensate between two 
parts of an optical lattice separated by an energy barrier is theoretically investigated. We show 
that by a pulsewise change of the barrier height, it is possible to switch between tunneling 
regime and a self-trapped state of the condensate. This property of the system is explained 
by effectively reducing the dynamics to the nonlinear problem of a particle moving in a double 
square well potential. The analysis is made for both attractive and repulsive interatomic forces, 
and it highlights the experimental relevance of our findings. 



The term Macroscopic Quantum Tunneling (MQT) usually refers to a property of a macro- 
scopic system, or part of it, that performs tunneling through a potential barrier like a single 
particle. After the first experimental discovery of this phenomenon [1] as a tunneling escape 
from a mctastable state in Josephson Junctions (JJ), there have been a number of similar 
realizations in completely different physical systems, such as liquid Helium [2] and nanomag- 
nets [3]. Frequently MQT is associated with the tunneling between different states of the 
system with no reference to a spatial energy barrier [4-7] , but there are also several examples 
where macroscopic objects, such as vortices [8,9] or fluxons [10] in JJ's, or magnetic domain 
walls [11], tunnel through a potential barrier. 

Recently [12,13], following earlier theoretical predictions [14], it has been found that a 
BEC trapped in an harmonic well potential of mesoscopic length behaves like a single JJ: 
for nonzero initial imbalance of the number of atoms in different wells, Josephson oscillations 
are present in the system, i.e. the condensate tunnels back and forth through the barrier. 
The only difference with respect to superconducting JJ's is that, for large initial imbalance, 
the condensate is mainly trapped in one of the wells, producing what is called Macroscopic 
Self-Trapping (MST). 
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Our aim in this Letter is to suggest the experimental realization of a BEC in an optical 
lattice, which is engineered in such a way to mimic two weakly coupled chains of JJ's (see 
Fig. Using such an experimental set-up, we demonstrate the feasibility of the efficient 
control of a switch between tunneling (MQT) and trapped (MST) states of the system. We 
show that our problem reduces to the Gross-Pitaevskii equation (GPE) [15] in a double square 
well, which displays very different properties from the previously considered double harmonic 
well potential [12, 14, 16-18]. Specifically, we show that for both attractive and repulsive 
nonlinearities the stationary solutions describing the MQT and MST regimes are characterized 
by very close energies in a wide range of the nonlinearity parameter. This property itself allows 
one to switch from the oscillatory tunneling regime to the trapped one and back via a simple 
pulselike adiabatic change of the energy barrier. Our results are broadly applicable and open 
the way to the experimental study of these phenomena in BEC dynamics. 




Fig. 1 - Schematics of the suggested experimental setup. The optical lattice is supplemented by large 
energy barriers from both sides and a small one in the middle (solid curve) . The condensate is initially 
loaded mainly into the right part of the optical lattice (the dashed line represents particle density). 
The inset shows the reduction of the problem to the particle motion in a double square well potential 
(details are given in the text). 

We start from the following, one-dimensional Hamiltonian of a BEC in an optical lattice: 

dip h 2 d 2 ip , . 2h 2 a s , „ 

where m is atomic mass, a s is the scattering length (a s < corresponds to attractive atom- 
atom interactions and a s > to repulsive interactions) and a± = ^/h/muj± is the transversal 
oscillation length, which implicitly takes into account the real three dimensionality of the 
system [19], loj_ being the transversal frequency of the trap. The optical lattice potential is 

V(x) — v cos 2 (k,Lx) for \ki,x\ > tt/2 
V(x) = (v + Vo) cos 2 (k L x) for \k L x\ < tt/2, (2) 

where is the wavenumber of the laser beams that create the optical lattice and Vq is the 
height of the additional spatial energy barrier placed in the middle of the optical lattice. 
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Besides that, Dirichlet boundary conditions with -0(±L) = are chosen in order to describe 
the large confining barriers at both ends of the BEC. These boundary conditions could be 
realized experimentally by an additional optical lattice with larger amplitude and larger lattice 
constant, as shown in FigQ 

Introducing a dimensionless length scale x = 2kLX and time t — Est/h, where Eb = 
8Er = Ah 2 k 2 L /m and Er is the recoil energy [20], we can rewrite JU as follows 



.9* _ ld 2 ^ 
% ~dt ~ "Ydx 2 



V(x)* + g\*\ 2 *, (3) 



where the normalized wave-function, J \^(x)\ 2 dx = 1, is introduced [21]. The dimensionless 
potential V still has the form J2J with the following dimensionless depths of the optical lattice 

v t7 Vq Na s 

v = w 9 = k^r (4) 




Fig. 2 - (color online) Numerical simulations of Eq. 10 : the left graph represents the transition from 
a self-trapped state to the macroscopic tunneling regime, while the right graph describes the inverse 
process. The insets in both graphs show the variation of the energy barrier necessary to realize the 
switching between the different regimes. 

As seen from the left panel of Fig. if one prepares the condensate in a self-trapped 
state it remains there until we apply the pulselike time variation of the barrier displayed in 
the inset. After that action, the condensate goes into the oscillating tunneling regime. On 
the other hand, preparing the condensate in the oscillating tunneling regime (right graph 
in Fig. [5J one can easily arrive at a self-trapping state by varying again the energy barrier 
in the middle as displayed in the inset. Let us mention that, as far as the energy of the 
barrier is changed adiabatically, the total energy of the condensate does not vary, i.e. the self- 
trapped and tunneling oscillatory regimes have the same energy. This is quite different from 



4 



EUROPHYSICS LETTERS 




Distance 



Fig. 3 - Stationary profiles described by expressions (I I 1 1 and il l 2t . The solid line represents symmetric 
and asymmetric solution profiles for attractive nonlinearities, while the dashed curves describe the case 
with repulsive interactions. The inset shows the relative energy difference between asymmetric and 
symmetric stationary solutions (solid and dashed lines correspond again to attractive and repulsive 
nonlinearities, respectively). 



what happens in a double harmonic well potential [12,14,16-18]. The point is that, in the 
double harmonic well, the asymmetric stationary solution is characterized by a smaller energy 
than the symmetric solution and this difference increases sharply with increasing nonlinearity. 
Hence, a drastic energy injection is required in order to realize the transition between the 
two regimes; whilst in our case the transition is simply achieved only by varying pulsewise 
the energy barrier. Below we argue that this happens because our case effectively reduces 
to the case of a double square well potential (see the inset of Fig. ^ an d the reduction 
procedure below) for which asymmetric and symmetric stationary solutions carry almost the 
same energies in a wide range of the nonlinearity parameter. 

Now we proceed to reducing Eq. © to a Discrete NonLinear Schrodinger equation (DNLS) . 
We discretize it via a tight-binding approximation [22-24], representing the wave function ^(x) 
as 

9(x) = Y6 j $ j {x), (5) 



where &j(x) is a normalized isolated wave function in an optical lattice in the fully linear case 
9 = and could be expressed in terms of Wannier functions (see, e.g., [25]). For clarity, we 
use here its approximation for a harmonic trap centered at the points rj — jir(\j\ + 1/2) 

varies from 1 to n, the number of wells). In the context of the evolution equation (J3J 
&j(x) has the form 

f° r \j\ 7^ 1 1 an d one should substitute v by v + Vq in the above expression in order to get an 
approximate formula for the wave function for \j\ = 1. 
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Assuming further that the overlap of the wave functions in neighboring sites is small, we 
get from © the following DNLS equation for the sites \j\ ^ 1 



dt 



-Q{<t>. 



j+i + (Pj-ij 



while for |j| = 1 we have 



(7) 



ih- 



3±1 



dt 



-Q0±2 - + Ui\4>±i\ 2 4>±i, 



(8) 



where we assume pinned boundary conditions. The constants Q, Q±, U and U\ are easily 
computed from the following expressions ^ 0): 



dx dx 
dx dx 



dx. 



+ {v + Vb)cos 2 (i/2)$i$_i 



dx, 



U = g J <f> 4 j dx ~ U! = g J <f> 4 ±1 dx 



(9) 



In order to characterize the solutions of Eqs. @ and |JS}, we follow the same procedure used 
in Ref. [26], which goes through a continuum approximation. Assuming that <f>\ — 4>-\ we 
finally arrive at 



Q dt 



dp 



w{j)m+Rm)?m, 



(10) 



where now j is a continuous variable, W(j) is a double square well potential with a barrier 
height w = 2(Q — Qi)/Q and width / = 1, (p(j) obeys pinned boundary conditions 4>{j = 
±L) = (2L is a width of a double square well potential) and the nonlinearity parameter is 
given by R — U/Q. 

Summarizing, we have reduced the initial problem, GPE with optical lattice and barrier 
potentials to a DNLS equation and then this latter again to a GPE with double square well 
potential. The reason for doing this, is to get rid of the optical lattice potential and to reduce 
our problem to the GPE with a double square well potential, for which one can easily find 
exact stationary solutions. They are sought as 4>(t,j) — <£(j) exp(— i{3t) with a real-valued 
function <&(x) found in terms of Jacobi elliptic functions [27], in the case of attractive atom- 
atom interactions R < 



-L < x < -I 
I < x < L 
-I < x < I 



$ = J 4cn[7 A (j + L)-K(k A ),k A ], 

$ = Bcn[7bG-L) + K(kB),k B ], 
$ = C dn[ 7c (j -j ),k c ], 



(11) 



with the parameters given in terms of the amplitudes by 
A 2 , , A 2 



Ya = P + 



R\ 



k 2 A 



lc = w -f 3 



2(A 2 + \R\py 
C 2 



2\R\ 



w~(3~C 2 /\R\ 
w - f3-C 2 /2\R\' 



B 2 



2{B 2 + \R\PY 
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while, in case of repulsive interactions R > 0, one has the stationary solution written in the 
form 



-L < x < -I 
I < x < L 
-I <x <l 



$ = A sn[ 7A (j +L),k A ], 
$ = B sn[ 7B (j -L),k B ], 

a> = C/cn[ 7c (j-j ),k c ], (12) 



where 



Q j4 q A. q B cy B 

lA=P- 77T75T> k A = 7^77, To, 7b = P ~ 7777*, fc B = 



2|i?|' 2|i?|/3-A 2 ' /jB ^ 2|i?|' * 2|i?|/3-5 2 ' 

7c = 2 f w -/3+SV fcc 



\ 2 w-/?+C 2 /2|fl| 

' c " w -f3 + cy\R\ ■ 

Here K denotes the complete elliptic integral of the first kind, and, by construction, the above 
expressions verify the vanishing boundary values in j = ±L. 

The solutions are then given in terms of five parameters (A, B, C, f3, jo), which are 
determined by the four continuity conditions in j = ±1 and the wave function normalization 
condition J dj& 2 (j) = 1. In both the repulsive and attractive case, one has symmetric and 
antisymmetric solutions, and the asymmetric solution appears above a given nonlinearity 
threshold value. In Fig. |2| we plot the profiles of the lowest energy symmetric (attractive case) 
and antisymmetric (repulsive case) solution, the asymmetric solutions for both the repulsive 
and the attractive cases, and the relative energy differences AE — 2(E a — E s )/(E a + E s ) 
between asymmetric E a and symmetric E s cases. As seen from the inset, these energies are 
very close in the nonlinearity range 0.018 < \g\ < 0.03 (note that the numerical simulations 
presented in Fig. |2]are made for g = —0.025) and hence it is easy to switch from the tunneling 
regime to the self-trapped state and back again. 

In order to get an idea about a possible experimental realization of these effects, we choose 
7 Li, which is characterized by an attractive atom-atom interaction [28]. With the realistic 
experimental parameters uj± = 2ir x 30Hz, a s = — 1.4nm, = 7.4 • /im -1 , from the formula 
for nonlinearity g — Na s /oj_ki, one gets that the total atom number needed to access values of 
\g\ around 0.025 is N w 10000. While in case of 23 Na [29] with repulsive forces (a s — 4.9nm), 
the atom number should be N w 1000. Increasing the number of wells (system size) n times, 
one should decrease the number of atoms n 2 times in order to observe the predicted effect. 

To summarize, we predict an effective switching scenario from an oscillatory tunneling to 
a self-trapped regime as a novel macroscopic quantum transport effect to be realized with 
ultracold atoms. We have shown that the problem effectively reduces to the particle motion 
in a double square well potential in contrast to earlier studies dealing with double harmonic 
wells. This difference guarantees the possibility of the switch from an oscillatory tunneling to 
a self-trapped state via a pulse-like change of the central potential barrier. We have derived 
typical ranges of physical parameters, for both attractive and repulsive interatomic forces, in 
order to suggest a ready-to-implcmcnt experimental verification. 
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